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FREQUENCY RESPONSE OF A VAPORIZATION PROCESS 
TO DISTORTED ACOUSTIC DISTURBANCES 
by Marcus F. Heidmann 
Lewis Research Center 

SUMMARY 

The open- loop frequency response properties of a drop vaporization process were 
numerically evaluated for acoustic disturbances with harmonically distorted wave shapes 
(wave shapes which deviated from pure sinusoidal forms). A single- drop model was 
used which allows for oscillations in drop mass and temperature. Calculations were 
made for n- heptane droplets vaporizing in oxygen \mder simulated unstable rocket com- 
bustor conditions. Response factors expressing the mass vaporized in phase and out of 
phase with the pressure oscillations were evaluated over a broad range of frequencies 
for a variety of harmonically distorted disturbances. 

For sinusoidal oscillations, typical response properties of the vaporization process 
were obtained with the single- drop model. The in-phase response factor increased from 
zero at low frequencies to a peak value at some intermediate frequency and then de- 
creased to a large negative value at high frequencies. Harmonic distortion was fotmd to 
substantially affect these in-phase response properties. In some instances the peak re- 
sponse was increased by an order of magnitude above that for sinusoidal disturbances. 

In general, the response factor increased by a constant increment for all frequencies at 
and above that for peak response. At low frequencies the response always decreased to 
zero. Effects of comparable magnitude were found for the out- of- phase response factor, 
but they occurred at different frequency conditions. 

The effect of harmonic distortion on the peak in-phase response factor was equiva- 
lent to that reported for a vaporization model which neglected the frequency- dependent 
behavior associated with drop mass and temperature oscillations. That study concluded 
that harmonic distortion can cause or contribute to the instability in rocket combustors. 
The present study for the frequency- dependent process gives the same conclusions. The 
present analysis also shows that harmonic distortion can substantially broaden the fre- 
quency range over which the vaporization process can support imstable combustion. 



IINTTRODUCTION 


A previous study (ref. 1) showed that the manner in which a vaporizing drop re- 
sponds to an acoustic disturbance is highly dependent on the wave shape of the disturb- 
ance. Harmonic distortion (deviations from sinusoidal wave shapes) was found to sub- 
stantially affect those response properties which control the stability of combustion 
systems. The study qualitatively showed that the harmonic distortion in a disturbance 
can cause or contribute to combustion instability problems in liquid- propellant rocket 
engines and provided some insight into the design of stable combustors. These results 
were obtained by neglecting the frequency- dependent behavior of the vaporization pro- 
cess. The present study extends the analysis into the frequency domain. The purpose 
of the study is (1) to verify the deductions of the previous study with a frequency- 
sensitive vaporization model and (2) to relate to combustor stability problems any 
changes in frequency- dependent behavior caused by wave distortions. 

The vaporization model used in the previous study of response properties (ref. 1) 
considered vaporization rate to be simply proportional to the drop Re 3 molds number to 
an exponential power. Such a quasi- steady model will provide no frequency- dependent 
behavior, whereas the response of the vaporization process has been shown to be highly 
frequency dependent (refs. 2 and 3). The sensitivity to frequency is caused by mass 
and thermal time lags; that is, a change in the vaporizing environment does not cause 
an instantaneous change in the drop mass and temperature. A model which includes a 
heat and mass balance for the vaporizing liquid drop is needed to give these time- lag 
properties. The present study uses such a model. 

The dynamic behavior of this vaporization process (model) is niunerically analyzed 
in this study. The specific process examined is n- heptane droplets vaporizing in gas- 
eous oxygen under simulated rocket combustion conditions. Calculations are made of 
the instantaneous vaporization rates caused by assumed acoustic disturbances (open- loop 
analysis). Evaluations are made over a broad range of frequencies for a variety of har- 
monically distorted distiirbances. The calculated oscillations in vaporization rate are 
expressed in terms of response properties which describe the relation between the 
vaporization rate and acoustic oscillations. The results are discussed within the con- 
tent of the previous studies and with regard to their application to rocket combustors. 

The response properties used in this study to relate vaporization rate oscillations 
to harmonically distorted acoustic oscillations require some clarification. These re- 
sponse properties are expressed as two quantities: the in-phase and out- of- phase re- 
sponse factors. These response factors are the in-phase and out- of- phase components 
of the vaporization rate oscillations relative to the acoustic pressime oscillations, where 
both oscillations are given as fractional perturbations about the mean value of the vari- 
able. The components are extracted directly from the distorted oscillations by a corre- 
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lation procedure which accounts for any contribution caused by harmonic content in the 
oscillations. Justification for such nonlinear evaluations of the response factors is 
given in reference 1. The basic significance of the response factors to combustion in- 
stability problems has evolved and been adopted from linear stability analyses. 

In linear (small sinusoidal) analyses the relation between two periodic variables is 
readily expressed either in terms of a phase angle and amplitude ratio or in terms of 
in-phase (real) and out-of-phase (imaginary) components. Both notations are interre- 
lated, but the in-phase and out- of- phase components have received most attention in 
combustor stability analyses. In stability analyses of combustor systems the boundary 
conditions characterizing the combustion process are specified by the in-phase and out- 
of-phase components of the burning rate oscillations relative to the pressure oscilla- 
tions. Both components (response factors) are needed to determine the stability char- 
acteristics of any specific combustion system. 

Systems analyses such as reference 4 show that the periodic behavior of any sys- 
tem is, in general, related to the properties of the response factors. The amplitude of 
sustained oscillations will increase with an increase in the in-phase response factor, 
and sustained oscillations will usually occur when the in-phase response factor exceeds 
a threshold value of about unity. The in-phase response factor, therefore, can be used 
as a qualitative index of the stability of a rocket combustor system. The out-of-phase 
response factor mainly controls the frequency of sustained oscillations. The frequency 
increases with an increase in the out- of- phase response factor. 


SYMBOLS 


c speed of sound, m/sec; ft/sec 

H* drop enthalpy, H/H 

J out-of-phase response factor, eq. (10) 

m* drop mass, m/w 

n harmonic order 

P* environmental pressure, P/P 

P* drop vapor pressure, P.^/P 

Pj^,P 2 harmonic coefficients for pressure distiirbance (pj^ < 1.0), dimensionless 
q* heat transfer rate, q/wH 

(R. in-phase response factor, eq. (9) 

r drop radius, iJ.m 
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initial drop radius, iim 
t time, sec 

u transverse gas velocity, m/sec; ft/sec 

Av steady relative axial drop velocity, m/sec; ft/sec 
w* drop vaporization rate, w/w 

y ratio of specific heats 

9 phase relation between velocity and pressure disturbance 

X* heat of vaporization, A/H 

T characteristic drop lifetime, m/w, sec 

cp phase relation between harmonic components of disturbance 

(jo frequency, rad/sec 

Superscripts: 

( )' primed quantities denote dimensionless perturbations, X’ =- ~ 2 ^ 

_ X 

( ) barred quantities denote mean values 

OPEN-LOOP MODEL 

The response properties of the vaporization process are obtained by an open- loop 
analysis by which the perturbations in vaporization rate caused by assmned acoustic 
oscillations are evaluated and correlated. 


Vaporization Process 

The vaporization model used in this study closely follows that developed in refer- 
ence 5 for vaporization- limited combustion in liquid-propellant rocket combustors. The 
model is applied in a manner similar to that used in reference 3 for a linear analysis 
of the frequency response properties of the vaporization process. The linear analysis 
showed that a single- drop model can adequately characterize the entire vaporization 
process with its multiplicity of drops. The single- drop model must, however, account 
for the continuous injection of propellant which occurs in combustors and which is an 
important source of frequency- dependent behavior. A steady liquid flow into a single 
drop simulates the dynamic effect of continuous propellant injection into a combustor. 
The single-drop model used previously and in this study is comparable to experimental 
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models where a single drop is suspended and continuously supplied by a h3T30dermic 
needle . 

The following sketch depicts the analytical model; 
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A vaporizing drop of radius r is continuously supplied by a steady flow rate w. The 
heat transfer rate q* to the drop causes a vaporization rate w* . The total heat and 
mass transferred establishes the instantaneous drop mass m* and the drop enthalpy 
H* . 

The analysis of response properties is restricted to equilibrium conditions of vapor- 
ization. For vaporization in a steady environment, an equilibrium condition of heat and 
mass transfer exists, where none of the vaporization variables change with time. At 
equilibrium there is a specific drop radius for each steady injected flow rate which sat- 
isfies a given heat and mass balance. A comparable condition exists for vaporization 
in a perturbed environment where the flow field varies periodically at a constant fre- 
quency and amplitude. In this case the drop radius and the heat and mass transfer vary 
periodically about some mean value. At equilibrium the oscillations about the mean do 
not vary from one period to the next. This condition of periodic equilibrium is specif- 
ically analyzed in this study. 

An analysis is made of n- heptane vaporizing in gaseous oxygen at simulated rocket 
combustor conditions. The vaporization equations and physical properties given in ref- 
erence 5 are used. The specific equations used in this study (expressed in normalized 
parameters and including physical properties) are developed in appendix A and given by 
the following relations: 
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(1) Vaporization rate 


w* 


(Av/100)^/^ 

^ (ryiOO)^/^ 


(P* 


1 + 




( 1 ) 


where Cj^ equals 60 and 33, respectively, for Av given in meters per second and feet 
per second 

(2) Heat transfer rate 


q* 


c (AVIOO)^^^ 

^ (ryiOO)^/^ 



( 2 ) 


where C 2 equals 363 and 200, respectively, for Av given in meters per second and 
feet per second 

(3) Enthalpy balance 


=— (q* - w*A*) (3) 

dt m* 


(4) Mass balance 


= 1 - w* (4) 

dt 

The vapor pressure and heat of vaporization properties of n- heptane were curve fitted to 
give 


p* (5) 

X* = 0. 65 - — (6) 

2 

The coefficients for the vaporization- rate and heat- transfer- rate expressions are given 
for a condition of equilibrium vaporization in a steady environment. The condition used 
is a drop radius r^ of 100 micrometers, a steady relative axial velocity between the 
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drop and gas Av of 30. 5 meters per second (100 ft/sec), and a combustor pressure 
of 2068 kN/m (300 psia). At these conditions w* , q* , and H* are equal to unity. 

The variable m* is actually the characteristic drop lifetime (m* = m/w = t) and is not 
equal to unity. 


Acoustic Oscillations 


The acoustic oscillations used to perturb the vaporization process are identical to 
those assumed for some specific evaluations of reference 1 with which the results of 
this study will be compared. The oscillations in pressure and acoustic particle veloc- 
ity which describe the acoustic disturbances are assumed to be similarly distorted. 
The distortion consists of harmonic components which are added to a fundamental os- 
cillation. Distortion is varied by changes in the magnitude and phase of the harmonic 
components relative to the fundamental. Two types of distorted oscillations are exam- 
ined. 

(1) Second harmonic distortion 


Po 

P* ' = Ph cos wt + — cos(2ojt - cp ) 
Pi 




u =- 


CPi 


Po 

(cos wt - 9) + — cos(2wt - cp - 0) 

Pi 


(7) 


(2) Multiharmonic distortion 


P* ' = 'y ^ Pj cos(ncot - cp) 
n=l 


> 

oo 

u = - 2^ Pj cos(ncijt - cp - 9) 

J 


( 8 ) 


The phase angles cp and 9 express, respectively, the phase relation between the har- 
monic components and the phase relation between velocity and pressure. Equation (8) 
characterizes high- amplitude oscillations observed in unstable rocket combustors, where 
the amplitude of the harmonic component usually decreases exponentially with harmonic 
order, that is, p^ = p^. 
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Response Factors 


Two response factors (more fully discussed in the INTRODUCTION and in ref. 1) 
are evaluated to express the mass vaporized in-phase and out- of- phase with the pres- 
sure oscillation. The correlation procedure used to evaluate the response factors is 
given by the following expressions: 

(1) In-phase response factor (R 


(R = 



P* 'w* ’ do)t 



(2) Out- of- phase response factor J 


/ 2n 



METHOD OF ANALYSIS 
Numerical Solutions 


(9) 


( 10 ) 


The response factors were evaluated by niimerical procedures for the finite ampli- 
tude disturbances. A desk-top calculator (Hewlett-Packard 9100A) with a capacity equi- 
valent to 140 program steps and six storage registers was used. Solutions for the re- 
sponse factors were obtained by a convergence procedure involving cyclic repetition. 

The flow diagram shown in figure 1 describes the procedure. 

The calculation procedure is such that an initial drop radius r^ and initial enthalpy 
Hq for a steady environment are known, but the mean drop radius r and mean enthalpy 
H for a periodic environment are not known until valid solutions for (Jl and J are ob- 
tained. The enthalpy differs negligibly for steady and unsteady conditions, but r can be 
considerably smaller than r^. The initial drop radius r^ is comparable to the drop 
size delivered by the propellant injector during combustion instability. Some of the re- 
sults are presented as a function of this pseudo- injected drop size. 
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Figure 1. - Flow diagram showing convergence prxedure. 


Linear Analytical Solutions 

Linear anal 5 rtical solutions for the response factors were evaluated for comparison 
with the computer- generated numerical solutions at small pressure amplitudes. Linear 
anal 5 dical solutions for the frequency response properties of n- heptane have previously 
been reported in reference 3. In that study correction terms for simultaneous heat and 
mass transfer of the drop surface and the variation of X (heat of vaporization) were ne- 
glected. Appendix B of this report presents a linear analysis which follows the method 
of reference 2 but includes these additional variables. 

The transfer function which applies specifically to this study as derived in appen- 
dix B is given by 
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w* ' ^ Q_ g22 2m* S /l - 0.249 m*s \ 

p*' ‘ 1 + 2m*S \1 + 0. 178 m*S/ 

The response factors (R and J are the real and imaginary parts of w* '/P* 
where S = icu. 


RESULTS AND DISCUSSION 

Harmonic distortion in the acoustic oscillations affected the frequency response pro- 
perties of the vaporization process in a similar manner for both types of distortion. 

This general effect of distortion on the frequency response is described first. Some re- 
sults follow for specific variations in second harmonic and multiharmonic distortion 
similar to those used in the quasi- steady analysis of reference 1 . 


General Effect of Distortion 

A general effect of harmonic distortion on the frequency response observed for all 
t 3 q>es of distortion is illustrated in figure 2. The in-phase response factor (fl is shown 
as a function of frequency for both sinusoidal oscillations (P 2 /Pj^ = 0) and for oscillations 
distorted by second harmonic content (P 2 /P 1 = 0*2) under otherwise identical conditions. 
The sinusoidal oscillations give the t 3 q)ical frequency response properties reported pre- 
viously (refs. 2 and 3). The response increases from zero at low frequencies to a peak 
value at intermediate frequencies and then decreases to a constant negative value at high 
frequencies. Although the response always decreases to zero at low frequencies, fig- 
ure 2 shows that harmonic distortion generally elevates the response curve above that 
for sinusoidal disturbances. Distortion increases the response factor by a relatively 
constant increment for frequencies at and above that for the peak value of the response 
factor, implying an additive rather than amplifying effect of distortion on response. The 
amoimt of the incremental increase is nearly identical to that obtained for comparable 
conditions of distortion in the quasi- steady analysis of reference 1. 

General agreement is exhibited between the results of this study and that of refer- 
ence 1 concerning the increase in the in-phase response factor caused by harmonic dis- 
tortion. The conclusions stated in reference 1 regarding the effect of distortion on the 
peak response appear valid and the consequences of changes in the in-phase response 
factor on combustor stability extensively discussed in reference 1 will not be reiterated 
in this report. Some additional conclusions regarding frequency- dependent behavior, 
however, can be drawn from the present study. In particular, figure 2 shows that har- 
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Frequency. Hz 


Figure 2. - Characteristic effect of wave distortion on in-phase response factor. Con- 
ditions: initial drop radius rQ, 100 |im; relative axial Mach number Av/c, 0.02; 
pressure amplitude 0.02; ratio of harmonic coefficients ^2!^' 
torted disturbances and 0 for sinusoidal disturbances. 

monic distortion can broaden the frequency range over which any vaporization process 
can cause instability. For example, if a value of (R greater than 0. 6 will cause insta- 
bility in a particular combustor, figure 2 shows that the instability can only occur at a 
frequency of 100 hertz with sinusoidal oscillations, whereas it can occur at any frequen- 
cy larger than 20 hertz with the distorted oscillations. 

Harmonic distortion also affects the out-of-phase response factor y, but the effect 
is less readily generalized than that for the in-phase response factor (R. The effect of 
distortion in ^ is best described by examining the response properties in the complex 
plane of (R and ^ shown in figure 3. The response curve for sinusoidal oscillations 
(P 2 /P 1 = 0) is t 3 q)ical of that used as the boundary condition for vaporization- limited 
combustion to obtain stability solutions for rocket combustor systems (ref. 4). The 
value of y is shown to vary between 0. 2 and -0. 85. The value of y is zero when (R 
is a maximum. For the condition of P 2 /Pj^ = 0. 2 (identical condition for that of fig, 2) 
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Figure 3. - Characteristic effect of wave distortion on frequency response properties in 
complex plane. Conditions: initial drop radius r^, 100 lam; relative axial Mach num- 
ber Av/c, 0,02; pressure amplitude 0.02; ratio of harmonic coefficients P 2 /P 1 , 
0.2 for distorted disturbances and 0 for sinusoidal disturbances. 


the distortion mostly affects y at frequencies less than that for the peak value of (R, 
and the effect is substantial. The most important effect of distortion on the response 
curve in the (R, ./ plane is the translation in the <R- coordinate caused by the additive 
effect of distortion on (R noted in figure 2. The change in however, is not insigni- 
ficant with regard to instability in rocket combustors. Any increase in the area encom- 
passed by the response curve in the (R, y plane increases the number of combustion 
chamber configurations which will give periodic solutions (instability) with a particular 
vaporization process. 

The specific effect of distortion on the response curve in the (R, ./ plane may differ 
from that shown in figure 3 for other conditions. Generally, distortion increases the 
area encompassed by the response curve and increases the prospect for vmstable com- 
bustion. 

Another generalization can be drawn from the results of the analysis which is not 
specifically related to harmonic distortion in the acoustic disturbances. Rather, it is 
a general behavior related to the amplitude of the disturbance. Figure 4 shows the gen- 
eral effect of pressure amplitude on the in-phase response curve for both distorted and 
sinusoidal disturbances. In addition to the change in the magnitude of the response fac- 
tor, figure 4(a) shows that the frequency at which the response attains a peak value in- 
creases with an increase in pressure amplitudes. This effect was also noted in refer- 
ence 2. The peak response frequency increases because the mean drop size, which 
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Distorted disturbances, P 2 /P 1 = 0.2 


Sinusoidal disturbances, P 2 /P 1 = 0 





1 . I 1 Llil 

10 ^ 10^ 11 

Frequency, Hz 
(a) Variation with frequency. 



1 10 1 
Dimensionless frequency, o/f 


(b) Correlation with dimensionless frequency. 

igure 4. ~ Effect of pressure amplitude on response properties. Conditions; initial 
drop radius rQ, 100 pm; relative axial Mach number Av/c, 0.02. 





characterizes that of an entire spray, decreases with an increase in the amplitude of 
the acoustic oscillations. The smaller drop size can respond more quickly to a change 
in the environment and thus shifts the peak response to a higher frequency region. 

The average change in drop size with pressure amplitude is shown in figure 5. 

This change in drop size with pressure amplitude is relatively independent of the amount 
of harmonic distortion when the amplitude is expressed (as in fig. 5) as a root-mean- 
square value. Also shown in figure 5 is the characteristic drop lifetime r, or its equi- 
valent m* , which is proportional to the cube of the drop radius. The characteristic 
drop lifetime is used in the dimensionless frequency wr, which is provided by linear 



Figure 5. - Effect of pressure amplitude on mean drop size and characteristic time. Conditions.- 
initial drop size Tq, 100 (jm; relative axial Mach number Av/c, 0.02. 


analyses (appendix B and ref. 3) as a correlation parameter for frequency response 
properties. The response properties of figure 4(a) are shown as a function of this 
dimensionless frequency in figure 4(b). The parameter correlates the results such that 
peak response for all conditions occurs at an oSf of about 1.5. 

Figure 3(b) also shows the linear solution for the in-phase response factor, which 
can be compared with the numerical solution for low- amplitude sinusoidal disturbances. 
Appendix C shows that the numerical solutions converge on the anal 3 dical solutions at 
extremely low amplitude and provide a degree of confidence in the numerical solution 
techniques. This is examined more thoroughly in appendix C. 

The shift in the frequency for peak response with pressure amplitude noted in fig- 
ure 3 implies that disturbances of a certain amplitude may trigger instability in a com- 
bustor by virtue of the frequency shift rather than by a change in the magnitude of the 
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peak response. For example, the natural resonant frequency of a combustor may be 
much larger than the frequency for peak response at small pressure amplitudes. An 
increase in pressure amplitude, however, could increase the peak response frequency 
until it matches the resonant frequency and thereby cause instability. Rocket combus- 
tors are often subjected to dynamic stability tests by introducing disturbances of in- 
creasing amplitude and harmonic distortion in the combustion chamber. A certain dis- 
turbance may cause instability not only because the distortion increased the in-phase 
response factor, but also because the amplitude of the disturbance shifted the peak re- 
sponse frequency. Figure 4 shows that a pressure amplitude of 0. 1 will increase the 
frequency for peak response by a factor of 6 above that for infinitesimal disturbances. 


Additional Effects of Harmonic Distortion 

In reference 1, second harmonic distortion produced a maximum in-phase response 
factor for particular values of second harmonic content and pressure amplitude. The 
frequency response properties exhibit a similar behavior. Figures 6(a) and (b) show 
the effect of pressure amplitude on the in-phase response factor (R for disturbances 
distorted by a constant amount of second harmonic content (P 2 /P 1 = constant). Figure 7 
shows the effect of the amount of second harmonic content at a constant pressure ampli- 
tude. Maximum response is shown to occur at a pressure amplitude of about 0. 02 with 
a second harmonic content P 2 /Pj of 0- 8. This maximum response is about an order 
of magnitude larger than that for sinusoidal disturbances (P 2 /PX = 0)* 

The effect of a change in the harmonic phase angle <p and the velocity- pressure 
phase angle 6 is shown in the (R, > plane in figure 8. An increase in these phase 
angles usually decreases both (R and ./ and thus decreases the probability of incurring 
instability. 

Multiharmonic distortion (eq. (8)) affects the frequency response properties, as 
shown in figure 9 and 10. Harmonic distortion increases with pressure amplitude and 
affects the pressure wave shape, as shown by the insert in figure 9. As inferred in 
reference 1, the peak in-phase response factor is shown to increase with an increase 
in pressure amplitude and to approach assnnptotically a maximum value of about 1. 4. 

The peak response frequency also increases with pressure amplitude. The effect of 
pressure amplitude on the response properties in the (R, > plane is shown in figure 10. 
The general effect of multiharmonic distortion on these response properties is similar 
to that caused by second harmonic distortion. 
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In-phase response factor, 
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gure 8. - Effect of harmonic phase angle and velocity- pres sure phase 
angle on response properties. Conditions: initial drop radius r^, 
100 urn; relative axial Mach number Av/c, 0.02: ratio of harmonic 


In-phase response factor, . 






Figure 10. - Comparison of response properties for multiharmonic and 
sinusoidal disturbances. Conditions: initial drop radius r^, 100 pm; 
relative axial Mach number Av/c, 0.02; pressure amplitude 
0 . 1 . 


SUMMARY OF RESULTS 

The results of this numerical analysis of the effect of wave distortion on the re- 
sponse properties of a frequency- dependent n- heptane vaporization process are sum- 
marized as follows : 

1. The amplifying effect of wave distortion on the response properties obtained by a 
previous analysis using a quasi- steady vaporization model were confirmed by the 
frequency- dependent model. The in-phase response factor exhibited substantial in- 
creases in peak response as predicted by the quasi- steady analysis. At frequencies 
higher than that for peak response, harmonic distortion increased the in-phase response 
factor by an increment equal to that for the peak response. At low frequencies the re- 
sponse always approached a value of zero. The frequency for peak response increased 
with the root-mean- square amplitude of the pressure oscillation. The increase was 
relatively independent of the type of harmonic distortion. The out-of-phase response 
factor was mostly affected at a frequency less than that for peak in-phase response. 

2. The results imply that the effects of harmonic distortion deduced from the quasi- 
steady analysis with regard to increasing the probability of spontaneous and triggered 
instability in rocket combustors are equally applicable for the frequency- dependent pro- 
cess. In addition to these previous effects concerned with the magnitude of the response 
factors, this study shows important changes in frequency characteristics. Distortion 
can substantially increase the range of frequencies in which the vaporization process can 
support combustion instability and can shift the frequency for peak response to match 
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the resonant frequency of a combustion chamber. In general, distortion increases the 
number of combustor configurations which can be imstable for a particular combustion 
process. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 13, 1972, 

113-31. 
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APPENDIX A 


VAPORIZATION EQUATIONS 

In this appendix equations for the vaporization process are developed for drops 
which are at or near their equilibrixim temperature. The vaporization model and physi- 
cal properties given in reference 5 are used for this purpose. The propellant combina- 
tion is n- heptane and gaseous oxygen. For convenience, an English system of units 
generally following that of reference 5 is used in developing the normalized equations. 


Physical Properties and Special Symbols 


The physical properties of the n- heptane - oxygen mixture surrounding the drop are 
assumed to be unaffected by changes in the drop temperature. Specifically, they are 
evaluated for a drop temperature of 845° R (equilibrium temperature for a chamber 
pressure of 300 psi) and a mean gas temperature in the drop boundary of 2000° R. The 
physical properties evaluated at these conditions and the dimensional symbols used in 
this appendix are given by the following notations: 



H 

k 

M 

M 


m 

V 


P,P 


Pr 

q 

R 


specific heat of heptane vapor, 0.9365 Btu/(lb)(°R) 
specific heat of heptane liquid, 0.70589 Btu/(lb)(°R) 
specific heat of gas mixture, 0.641 Btu/(lb)(°R) 

molecular diffusion coefficient of gas mixture, 7. 492x10"^ (300/P) in.^/sec 
drop enthalpy, Btu/lb 

thermal conductivity of gas mixture, 6. 673x10"'^ Btu/(in. )(sec)(°R) 

molecular weight of gas mixture, 46.3 Ib/mole 

molecular weight of heptane, 100 Ib/mole 

2 — 

total pressure, Ib/in. ; P = 300 psi 

7-[4670/(Cp) T,J 

vapor pressure of heptane, 300 e ^ psi 

Prandtl nxunber, (c ) ii/k = 2.912, dimensionless 

heat transfer rate, Btu/sec 

universal gas constant, 1545 (ft)(lbf)/(°R)(lb)(mole) 
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Re Resmolds number, 2rUp^/ju = 3. rPU, dimensionless 

r drop radius, pm 

Sc Schmidt number, ii/Dp^ = 1.0188, dimensionless 

T gas mixture temperature, 2000° R 

Tg combustion gas temperature, °R 

drop temperature, °R 
t time, sec 

U magnitude of relative drop velocity, ft/sec 
u transverse gas velocity, ft/sec 

Av steady relative axial drop velocity, ft/sec 

w vaporization rate, Ib/sec 

Z, z correction terms for heat transfer, eq. (A6) 

X heat of vaporization of heptane, 392 - T^/2 Btu/lb 

p viscosity of gas mixture, 3. 033x10"® lb/(in. )(sec)(°R) 
density of heptane liquid, 1.8646x10"^ Ib/in. ® 

Pm density of gas mixture, 3. 76x10" '^(P/300) Ib/in. ^ 
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the vaporization rate can be expressed as 


w = 0,32X10-= 

\ 100 / \ 100 / 




1/4 


In 


iP- P- 


(A4) 


Heat Transfer Rate 


The heat transfer rate is given by 


q = 2ffkr[2 + 0. 6(Pr)^/^(Re)^/^](Tg - T^Z 

where Z = z/(e^ - 1) and z = jw(Cp) (T^ - T^zj/q. For the stated properties equa- 
tion (A5) becomes ^ 

q = 1. 65X10” ^®r[2 + 1. 67xlO”2(rPU)^/^j(Tg - T^Z 


(A5) 


(A6) 


By treating the bracketed term as for vaporization rate and by assuming a constant 
temperature difference T - of 4000° R, the heat transfer rate expression becomes 


q = 1.95x10 


-3 


Tr \l/2 / V \3/2 / p\l/2 


Vioo/ \100 


1 + 


iS 


1/4 


(A7) 


Combining the vaporization rate and heat transfer rate expressions to evaluate z 


gives 


z = 6. 15 hi I 


i P - P- 


(A8) 


The heat transfer rate expression can be expressed as 


q = 12x10' 


,-3 /av^'^^ / r i 

\iool \100/ \P/ 




ql/4 


In 


iP - P.. 


6.15 


(A9) 



Mass Balance 


The mass balance for a vaporizing drop being replenished by a constant flow rate 
w is given by 


— = w - w (AlO) 

dt 


Heat Balance 

The heat balance for a drop where the difference between the sensible heat carried 
by the vaporizing liquid and by the replenishing liquid is neglected is given by 

dt m 


Normalized Equations 


Normalized variables for the vaporization expressions are given by 


i w 
w* 

w 


m 


* _ 


P* = . 
V 


m 


j* =-3^ 
wH 


P* 

P 


(A12) 


X* =- 


w 


H 


H* =-5 
H 


J 


where w is the vaporization rate of an initial drop size r^ for steady conditions and 
Av = 100 ft/sec as given by 


v3/2 


w = 0.189x10" 


. 100 ) 


and m is related to the drop radius by 

4 3 

m = — Trr p, 
3 ^ 


(A13) 


(A14) 
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With these notations, the vaporization process at a mean chamber pressure of 
300 psi is given by the following expressions: 

^ 33 (Av/100)^ ^ ^ )l/2 

(ryiOO)^/^ 


q* = 200 (^v/100)^'^^ . 

(ryiOO)^/^ 


= A (q* - w* X* ) 
dt m 

P* ^e'^A7.8/H*) 
X* = 0. 65 - H* 



(A15) 

(A16) 

(A17) 

(A18) 

(A19) 
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APPENDIX B 


LINEAR ANALYSIS 

In this appendix a transfer function relating perturbations in pressure and vaporiza- 
tion rate is derived by linear analysis. 

The normalized equations for the vaporization process given in appendix A are used. 
The linearized forms of these equations when gas velocity effects are neglected are 
given by 


w* ' = — P* ' + — m* ’ - /3(P* ’ - P* ’) 
2 2 


(Bl) 


q* » p* . + i m* ' - /3 fl - 


Gof 


e®“ - 1 , 


|(p*' - P*0 


(B2) 


m*m* 'S = -w* ' 


(B3) 


^ in* H* ’S = q* ’ - X* ’ - w* ' (B4) 

X* 


P* ’ = 7. 8 H* ' 


X* » H* /2 , 

0. 65 - (H*/2) 

where 



(B5) 

(B6) 


(B7) 


(B8) 


For steady vaporization at a pressure of 300 psia, /3 = 1. 37, a = 0. 599, H* =1.0, and 
X* = 0. 155. 
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Evaluating the linearized equations at this condition gives 


w* ’ = -0.87 P* ' P*’ 

q* ' = 4. 2 P* ’ + m* ' - 3. 7 P* ’ 

m*m* 'S = -w* ' 

6. 45 m* H* 'S = q* ’ - - w* ’ 

P* * = 7. 8 H* ’ 

X* ’ = -3.23 H* ’ 

Combining these equations gives the transfer function 

w** ^ Q g22 2 m*S / l - 0.249 m*S \ 

P*' ' 1 + 2 m*S\l + 0. 178 m*S/ 


(B9) 

(BIO) 

(BID 

(B12) 

(B13) 

(B14) 

(B15) 
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APPENDIX C 


COMPARISON OF VAPORIZATION MODELS 

The frequency response properties of the vaporization process reported previously 
(refs. 2 and 3) were restricted to sinusoidal acoustic oscillations. A comparison with 
these previous results is made here (1) to bring the present analysis into perspective 
with the previous studies and (2) to establish a degree of confidence in the specific model 
and method of evaluation used in this study to examine the effects of wave distortion on 
response properties. 

The single- drop model used in this study attempts to simulate the entire vaporiza- 
tion process in a combustor with its multiplicity of drops and continuous injection. The 
response of the entire vaporization process to sinusoidal disturbances was evaluated for 
n- heptane drops in reference 2. Individual heating transients, accelerations, and vapor- 
ization histories for each drop within the process were evaluated. Some of the results 
of that study are shown in figure 11. The average value of the in-phase response factor 
(R is shown as a function of frequency for variations in drop size, final gas velocity, and 
pressure amplitude. For each condition there is a frequency for peak response. This 
frequency increases with a decrease in drop size and with an increase in final gas veloc- 
ity and pressure amplitude. The peak values of the response factor curves are rela- 
tively constant at a value of about 0. 7. 

Figure 12 shows comparable results obtained for the model and method of analysis 
used in the present study. (An increase in final gas velocity is comparable to some pro- 
portional increase in relative axial Mach number Av/c. ) The peak response, the fre- 
quency at which this peak occurs, and the effect of the variables on this peak response 
frequency are in relatively good agreement with those of reference 2. The model 
appears to adequately characterize the frequency response properties for the purpose of 
examining the effect of wave distortion on the response properties. 

In reference 3, analytical solutions for the frequency response properties of the 
vaporization process were derived by linear analysis. Such anal 3 dical solutions provide 
exact values of the response factors for infinitesimal sinusoidal oscillations. The solu- 
tions also provide the correlating parameter <ut (characteristic frequency) for the re- 
sponse properties. The parameter t is a characteristic drop lifetime defined by m/w, 
which for the numerical evaluations of this study is equal to m* . 

Figure 13 shows the in-phase response factor (R as a function of the characteristic 
frequency o)t obtained both by the linear analysis developed in appendix B and by 
numerical evaluations. Also shown in figure 13 (lower curve) is the result from the 
linear analysis of reference 3, where a correction term for simultaneous heat and mass 
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Response factor, dimensionless 



Drop 


radius, 


•"O' 


pm 

o 

50 

□ 

150 

A 

500 


(a) Effect of initial drop radius. Drop conditions: 
veixity, 30. 5 m/sec (100 ft/ sec); temperature, 
361 K (650° R). Gas conditions: pressure, 

2068 kN/m^ (300 psi); final velocity, 244 m/sec 
(800 ft/sec); oscillation peak-to-peak amplitude, 
0.2 of average gas pressure. 


Final gas 
velocity, 
m/sec 

□ 122 
O 244 
A 732 



(b) Effect of final gas veixity. Drop conditions: 
radius, 50 veixity, 30. 5 m/ sec (100 ft/ 
sec); temperature, 361 K (650° R). Gas condi- 
tions: pressure, 2068 kN/m^ (300 psi); oscilla- 
tion peak-to-peak amplitude, 0.2 of average 
gas pressure. 


Pressure 

amplitude, 

ratio, 

2pr 

□ 0.1 



(c) Effect of peak-to-peak pressure amplitude ratio. 
Drop conditions: radius, 50 pm; veixity, 

30.5 m/sec (100 ft/sec); temperature, 361 K 
(650° R). Gas conditions: pressure, 2068 kN/m^ 
(300 psi); final veixity, 244 m/ sec (800 ft/ sec). 



Figure 11. - Frequency response properties from reference 2 obtained with a multiple-drop transient vaporization model. 


N3 

CD 





Dimensionless frequency, cjt 
F igure 13. - Comparison of analyses. 


transfer and the variation in the heat of vaporization with drop temperatures were ne- 
glected. The present linear analysis includes these variables. 

A comparison of the two linear analyses shows that the additional variables increase 
the peak response to a value of about 0. 6 and increase the characteristic frequency for 
peak besponse to a value of about 1.3. These changes alleviate some of the concern ex- 
pressed in reference 3 about the magnitude and position of the peak response from that 
linear analysis. 

A comparison of the numerical solutions and the modified linear analysis in figure 3 
shows relatively good agreement of the response properties. As pressure amplitude de- 
creases, the numerical response properties converge toward the linear solution. The 
result gives added confidence in this study to evaluate the response properties. 
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